Laser beam filament at ion in fractal aggregates 
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We investigate filamentation of a cw laser beam in soft matter such as colloidal suspensions and 
fractal gels. The process, driven by electrostriction, is strongly affected by material properties, 
which are taken into account via the static structure factor, and have impact on the statistics of the 
light filaments. 



Light propagation in soft-matter such as colloidal sys- 
tems can be affected via a feedback mechanism by light- 
induced structural changes. In spite of the low powers 
needed to observe such phenomena [l[, the optical non- 
linear response of soft matter was generally overlooked 
with the exception of metal colloids |2| and liquid crys- 
tals 3. 

Here we consider the case of a colloidal suspension 
of dielectric particles dispersed in a solvent, character- 
ized by tunable interactions that can be responsible for 
long range correlations and even phase transitions such 
as gelation. The role of the structure of these materials 
in nonlinear optical processes is essentially unexplored. 
Large scale ordered structures can affect orientational Q , 
electroctrictive [JQ, and thermoforetic 5] mechanisms. 
For isotropic particles and negligible thermal gradients 
due to light absorption, the leading mechanism is ex- 
pected to be electrostriction: the particles are subject to 
forces induced by light intensity gradients, thus moving 
in the region with higher or lower intensity depending on 
the difference between their refractive index n s and that 
of the host medium |4|. In both cases, the optical 
beam experiences self-focusing, a process which has been 
described by a local Kerr law or intensity-dependent re- 
fractive index variation An = n^I, where is the Kerr 
coefficient and / the optical intensity [l[ . Recently, it has 
been shown that the nonlinear response of such materi- 
als is mediated by the static structure factor S(q) of the 
material 7] , thus turning out to be strongly nonlocal Q . 
While the strength of the nonlinearity (112) is given by 
the material compressibility, S(q) is strongly affected by 
the whole structure of the soft-material phase, e.g. by 
the presence of fractal aggregates which has impact on 
the nonlinear susceptibilities. 

In this Letter, we investigate laser beam filamenta- 
tion, a process which is well known to occur due to 
spatial modulational instability (MI) MI is ubiq- 

uitous in nonlinear science, occurring in fluidodynamics 
(Benjamin-Feir instability matter wave s 1 1(1 . plasma 
physic s Il2l . coherent oscillations in lattices 11], and spin 
waves [13 • In the context of optics, filamentation has 
been recently shown to be affected by the coherence prop- 
erty of light [l4| , by the artificial periodicity of nonlinear 



media ? an d by nonlocality ^| . Specifically, filamen- 
tation has been investigated in liquid crystals where the 
reorientational nonlinearity allows to tailor the nonlocal- 
ity pj| , and is at the origin of complex light dynamics 
|l8j | . It is therefore interesting to investigate how the 
dominant nonlinear mechanisms in soft matter, and par- 
ticularly in complex media such as fractal aggregates, 
affect the filamentation process. In addition, the con- 
sidered light induced perturbations can be reversible or 
irreversible (depending on the specific material relaxation 
dynamics); in the latter case the filaments can be used 
for the all-optical imprinting of optics circuitry. 

We start from the following model for the electrostric- 
tive response of soft matter p 
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where E is the complex envelope of the electric field 
at angular frequency uj, k = kofio = (uj/c)no is the 
wavenumber, no is the average refractive index, and 
= <9 2 + dy. Here p is the optically induced perturba- 
tion to the bulk particle density po, which is related to 
the optical intensity through the static structure factor 
p(q) = 'y e ZoS(q)I(q)/2nokBT, where tilde denotes 
Fourier transform, j e = po(de/dp) Po the electrostrictive 
coefficient 4] Zo is the vacuum impedance, ks is the 
Boltzmann constant, and T is the system temperature. 

MI is the process that triggers the generation of mul- 
tiple filaments due to exponential amplification of plane 
wave perturbations with small transverse wavenumbers 
Q 2 = Qx + Qy (changes in longitudinal wavenumber q z 
of the perturbation is commonly assumed to be negligi- 
ble due to the paraxial conditions), at the expenses of a 
pump beam E = Eo ex.p(i/3z). A linear stability analy- 
sis allows us to find the growth rate or MI gain 
G(q x ,q y ) as 
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I P S(q) 
P r Sq 



(qw y 



(2) 



where wq is a reference waist (of intensity), Zd = kw^ is 
the associated diffraction (Rayleigh) length, P = ttWqIq 
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FIG. 1: (Color online) Properties of HS colloids: (a) S(q) vs. 
normalized wavenumber qr s for different values of packing 
fraction n; (b) So ( nonlinear ity) and S2/S0 (nonlocality) vs. 
77. 
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FIG. 2: (Color online) (a) S(q) in the presence of aggregra- 
tion, for D = 2.2 and f = 10 4 r s (thick line), $ = 5x 10 3 r s 
(dashed line) and £ = 10 3 r s (thin line); (b) MI gain when 
P 10 3 P r and £ = w = 10 4 r s , with D 2.2 (thin line) and 
D — 1.2 (thick line). 



is the beam power, with Iq = uoEq/ '(2Zq) the peak in- 
tensity, and P r = irksT poe^cno/ (So^k^) is a reference 
power (P r = 3/nW for typical numbers r s = 10 nm, 
wo = 100 /mi, n s = 1.5, rih = 1.3, and a particle den- 
sity po = 10 21 m -3 ). The Kerr limit is recovered in 
Eq. © for a constant S(q) = So, where So = 5(0) is 
the ratio between the material and the ideal gas com- 
pressibility, and is proportional to the Kerr coefficient 
ri2 = Z 7 2 5o/4eon 2 ) /c j BTpo- Let us consider, first, the 
case of colloids made by hard spheres (HS) with radius 
r s and dielectric constant e s , for which the electrostric- 
tive coefficient is j e = 47rr 3 } e/ l (e s — e^)/(e s -h2e^), being 
the dielectric constant of the host medium (e.g. water). 
In this case, S(q), from 0, is displayed in Fig. ^1 for 
different values of packing fraction 77 of the spheres. At 
very low packing fractions the medium behaves as a Kerr 
medium with an approximately constant S(q). However, 
the curvature of S(q) at low wave- numbers increases with 
77, until a clear peak appears for tightly packed spheres. 
This means that material response becomes more nonlo- 
cal as spheres are more densely packed. At the same time, 
however, the nonlinearity (compressibility) decreases. A 
parabolic approximation S (q) = So+S2(qr s ) 2 — So+Kq 2 
allows us to show explicitly the dependence of nonlin- 
earity (So) and nonlocality (S2/S0) on 77 (see Fig. QJ)). 
However, we argue that nonlocal behavior can be hardly 
seen for HS. In fact, by defining a degree of nonlocality 
a 2 = K/(S w 2 ) = (S 2 /So)(r s /wo)\ for typical num- 
ber, e.g. wo = 100 fim and particle size r s = 10 nm, we 
find a 2 9* IO" 9 [S2/S0 = 0.1, see Fig. [TJ3]. Hence for HS 
a local model works quite well [l(, even at high packing 
fractions, as expected on the basis of the extremely short 
interaction range involved. 

The situation becomes drastically different when long 
range correlation sets in due to induced particle aggre- 
gation, as for instance in gelly or glassy arrested states. 
Either kinetic or structural arrest results in a peak in 
S(q) at q = due to enhanced compressibility and long 
range spatial correlations. One possible way to realize 
kinetic, out of equilibrium, gelation process is diffusion 
limited cluster aggregation (DLCA) induced by a strong 
and very short-range attraction between particles in ad- 
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FIG. 3: (Color online) Dependence of MI features on aggre- 
gate dimension £: peak gain Gza (a) and maximally amplified 
wavenumber (blue line, left ordinates) q m ax and bandwidth 
(red line, right ordinates) Ag (b) vs. £/r s . Here P = 10 3 P r , 
D = 1.5, iuo = 10 3 r s . (c-d) shows the same quantities as 
a function of fractal dimension D for fixed £ = lOi^o (wo = 
10 4 r s ), and powers Pi = 10 3 P r ,P 2 = 10 4 P r ,P 3 = 10 5 P r . 

dition to HS repulsion [2(j. This means that once two 
particles get close they stick irreversibily resulting into a 
so called fractal gel, i. e. a space filling network of inter- 
connected fractal clusters, characterized by well known 
expressions of the structure factor (see e.g. Eq. 4 in 
[2l|). Such S(q), parametrized by the long range correla- 
tion distance (aggregate dimension) £ and by the fractal 
dimension of the aggregates D, is reported in Fig. [2Ji. 
The corresponding gain profile [from Eq. (0)] exhibits a 
single band as shown in Fig. ^p. The fractal dimen- 
sion D and the aggregates characteristic length £ have a 
strong impact on the characteristics of MI. First, com- 
paring Figs. [2t and^i, it is clear that S(q) varies on a 
q— scale much smaller and the nonlocality is much higher 
in the presence of aggregates. 

Specifically, from the expansion S = So — Kq 2 we find 
a 2 = K/(Sowl) = D(D + l)(£/w ) 2 /6, showing that 
(for a fixed fractal dimension D) the nonlocality grows 
with the relative (to beam waist) cluster dimension (,/wq. 
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FIG. 4: (Color online) (a) Waist (FWHM) of the solitary 
wave solutions (solid lines) and associated index perturbation 
(dashed lines) vs. their power for D = 1.3 (thick) and D = 2.3 
(thin). Here wo — 10 4 r s and £ = 10u>o- (b) Intensity and 
index profile sampled at P = 70P r (D = 1.3). 



Conversely if the beam waist contains many aggregates 
instead of a spanning network, the situation becomes sim- 
ilar to that of many independent particles in the beam, 
and hence the system behaves locally. 

The dependence of the MI feature on the aggregate 
dimension are summarized in Fig. 01 for a typical value 
of power P = 10 3 P r (no qualitative change are observed 
by varying P). Importantly, in spite of the fact that ri2 
increases with £ (So grows due to the higher compressibil- 
ity), the enhanced nonlocality acts as a filter that frus- 
trates MI, and cause the peak gain C7, the MI bandwidth 
Aq or the maximally amplified spatial wavenumber q m ax 
to decrease dramatically with £. This is relevant when 
the colloidal sample undergoes phase transitions such as 
gelation, because the filamentation pattern is expected 
to change as the long range correlation associated with 
material structural changes developes. The dependence 
of the MI gain on the fractal dimension D is summarized 
in Fig. 03, d. Notably for D < 2 the MI process appears 
more pronounced and qmax strongly varies with D. We 
emphasize, however, that smaller D implies smaller So, 
and hence larger powers P r . 

Beyond the early stage of exponential amplification, 
MI is expected to give rise to filaments corresponding to 
solitary waves of Eq. (1) In this sense, the nonlocal- 
ity plays a key role, acting as a spatial filter that tends to 
flatten the induced index perturbation, thus stabilizing 
the single filaments. This is shown in Fig. ^ where we 
report the FWHM (full- width-half-maximum, henceforth 
waist) of the intensity and associated index change pro- 
files against filament power (the profiles are shown in Fig. 
Ht>). Unlike local Kerr media, the index perturbation as- 
sociated with the solitary waves turns out to be broader 
than the intensity distribution, and its waist tends to sat- 
urate. Once again, the fractal dimension plays a crucial 
role determining the size of filaments. 

In order to provide more direct evidence that the out- 
put filament pattern is strongly affected by the fractal 
dimension D, we resorted to the numerical integration 
of Eqs. (PJ. The key point is to recognize that in soft 
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FIG. 5: (Color online) Filament pattern for two different frac- 
tal dimensions: (a) D = 1.3 (distance z = OAzr); (d) D — 2.3 
(distance z — zr); (b,e) Corresponding index perturbations; 
(c,f) Corresponding Su g ht (squares) averaged over 40 noise 
realization; Sl — Su g ht/lF (solid line) and MI gain profile 
(dashed line) (P 1000P r , v v 10" 4 , w 10 4 r s , f w ). 

matter, the leading mechanism that triggers laser beam 
filamentation, is the presence of material density fluctu- 
ations. These are known to be ruled by S(g), and, in the 
spectral domain, can be written as p = S(q)f, where / 
is a 3D stochastic variable (white noise). To model this 
effect we generalize Eq. JIJ by resorting to the following 
stochastic PDE model 

dct 1 

i-^ + g + ^ = ; ^ = R(q 8 , q t )(\a\ 2 + D), (3) 

where we have introduced the dimensionless quantities 
R = S/S , (s,t) = (x,y)/w , C = z/zd, and v is 
a spectral noise term such that (v(s,tX)v(s' ,t f , CO) = 
Vp8(s — s')5(t — t')5(( — CO 7 the- brackets denoting sta- 
tistical average. The density is given by p — pN^, with 
Pn = [koZd(dn/ dp)]~ l , and ip is the sum of a nonlinear 
term and of a colored (through R = S/ So) noise. Equa- 
tions © are solved by the Heun method [23 , after a FFT 
pseudospectral discretization along s and t. We show re- 
sults obtained with a plane wave input a(( = 0) = ao 
(oq = P/AP r due to simple scaling arguments). 

At the early stages of the evolution the input plane 
wave develops a spatial modulation, which spectrally cor- 
responds to the MI gain bandwidth. At longer propa- 
gation distances narrow filaments are formed as evident 
by the snapshots (for a given noise realization) of the 
near-field intensity pattern \a\ 2 in Fig. I n the far- 
field such filaments are distributed along a circle due to 
the axial symmetry of the problem ("conical emission" 
[22| )• When the field propagate further, the filaments 
interact and form clusters (not shown) [l^- The fractal 
dimension D not only affects the MI process, but also 
the number of generated filaments. The global charac- 
terization of the filamentation process requires a statis- 
tical description; we introduce the Fourier transform of 
the autocorrelation of \a(x,y)\ 2 — (\a(x,y)\ 2 ) (the second 
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term allows to get rid of plane- wave background), aver- 
aged over a large number of noise realizations and along 
curves q = constant (due to radial symmetry). This 
quantity, which is the equivalent for light of the material 
S(q) and will be denoted as Su g ht(q), is shown in Fig. 
IS1 The peak positions of Su g ht are determined by the 
average filament distance. Additionally, so long as the 
filament positions and shape can be taken as indepen- 
dent we can write Su g ht — Sl(q)^f(q), where If(q) is 
a form factor, which is approximately the average line- 
shape of the filaments, and Sl{q) only depends on their 
statistical distribution. These quantities are shown in 
Fig. [5J assuming for If(q) a Gaussian profile (its width 
is the average filament waist). The first peak in Su g ht, 
or Sl, corresponds to the maximally amplified spatial 
harmonics of MI theory. 

In conclusion we have developed a theoretical descrip- 
tion of laser beam filamentation in fractal aggregates and 
soft matter in general. The long-range statistical proper- 
ties of multiple filaments reflect the material properties, 
and are strongly dependent by the fractal dimension of 
the aggregates. 
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